Electrostatics of Lipid Bilayer Bending 
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The electrostatic contribution to spontaneous membrane curvature is calculated within Poisson- 
Boltzmann theory under a variety of assumptions and emphasizing parameters in the physiological 
range. Asymmetric surface charges, either fixed with respect to bilayer midplane area, or with 
respect to the lipid-water area both induce curvature but of opposite sign. Unequal screening layers 
on the two sides of a vesicle (e.g. multivalent cationic proteins on one side and monovalent salt on the 
other) also induces bending. For reasonable parameters, tubules formed by electrostatically induced 
bending can have radii in the 50-100nm range, often seen in many intracellular organelles. Thus 
membrane associated proteins may induce curvature and subsequent budding, without themselves 
being intrinsically curved. Furthermore, we derive the previously unexplored effects of respecting 
the strict conservation of charge within the interior of a vesicle. The electrostatic component to the 
bending modulus is small under most of our conditions, and is left as an experimental parameter. 
The large parameter space of conditions is surveyed in an array of graphs. 

Keywords: membrane curvature, bilayer electrostatics, vesicles, Poisson-Boltzmann equation. 



INTRODUCTION 

The membranes bounding intracellular organelles are dy- 
namic structures, whose morphology the cell can reg- 
ulate. This happens most dramatically during mito- 
sis, when the nuclear membrane disintegrates, probably 
forming vesicles (Alberts, 1994). Proteins destined for 
secretion or targeted to the plasma membrane pass from 
the endoplasmic reticulum (ER) to the Golgi body and 
thence to their target; at each step vesicles or tubular pro- 
cesses with diameters in the 50-100nm range are involved 
in the sorting and transport (Rothman, 1994; Schekman 
and Orci, 1996). Similar remarks apply to endocytosis 
and endosomal sorting and maturation (Gruenberg and 
Maxfield, 1995; Trowbridge, et al., 1993). The ER cis- 
terna itself has a large tubular component (Terasaki et 
al., 1986) and there are tubular connections within the 
Golgi body (Rambourg and Clermont, 1990). Various 
treatments can enhance the tribulation of the membranes 
of Golgi and other organelles (Lippincott-Schwartz et al. , 
1991; Cluett et al, 1993). 

Numerous plausible mechanisms have been proposed 
for vesicle budding or tubule formation, but often there 
is insufficient physical detail for a quantitative assess- 
ment of their validity. For instance, are the building 
blocks of the clathrin cages intrinsically curved, or do ki- 
netic processes during budding favor five membered over 
six membered rings and thus effect closure (Shraiman, 
1996). The phase separation of wedge shaped lipids or 
membrane resident proteins, is also argued to play a role 
in vesiculation (Schekman and Orci, 1996), as is line ten- 
sion between the two phases (Lipowsky, 1993). Chemical 
modifications to the lipids on one side of a bilayer, such 
as phosphorylation of inositol lipids (de Camilli et al., 
1996), or cleavage of the acyl chains by phospholipases 
(Brown, 1996), are known to occur biologically and will 
promote membrane curvature. Less drastic modifications 
such as changes in solution pH or ionic strengths can have 



similar effects because many lipids are zwitterionic and 
their effective charge varies with solution conditions. Fur- 
thermore, since parameters such as pH, ionic strength, 
temperature, phospholipid pK a are not independent, it 
is impossible to separate and quantify all these effects. 

It has long been recognized that when surface charges 
are unequal, electrostatics will induce membrane curva- 
ture. In this Article we quantify a number of less obvious 
influences of electrostatics on membrane bending. Sym- 
metrically charged membranes (i.e., same surface charge 
on the two leaflets) will bend in response to asymmetri- 
cal screening. Thus membrane associated proteins such 
as the adaptins which mediate between the clathrin and 
the bilayer could cause bending, even if the proteins are 
globular while free in solution. Various plausible choices 
for the neutral surface with respect to which the surface 
charge is conserved during membrane bending can alter 
even the sign of the preferred curvature resulting from a 
given charge or screening asymmetry. Uncharged mem- 
brane lipid components can define the neutral surface and 
hence modulate the electrostatic response. Finally, if the 
interior of a vesicle is truly electrically isolated from the 
exterior, a small internal charge scaling with the area can 
define the spontaneous curvature. The simple descrip- 
tion of electrostatics we employ, the Poisson-Boltzmann 
equation, has been shown to give quantitatively accurate 
results for the arguably more complex problem of the 
binding energy between a charged protein and a bilayer 
(Tal et al. 1996). 

Even at the level of Poisson-Boltzmann theory, the 
electrostatic contribution to spontaneous curvature par- 
ticularly with asymmetrical (different on the two sides) 
and multivalent electrolytes have not been fully quanti- 
fied in the literature, and are potentially significant for 
the concentrated electrolytes encountered in the cell. Our 
interest was kindled by the ubiquity of tubules within the 
diameter range noted above, and the possibility of get- 
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ting such sizes by simple electrostatics. More attention 
has been paid in the physical literature to the electro- 
static contribution to the curvature energy (we will gen- 
erally take the total bending modulus as an experimen- 
tally determined parameter) rather than the spontaneous 
curvature. Low ionic strength (which enhances the elec- 
trostatic contribution to the bending modulus) was em- 
phasized when the spontaneous curvature was computed, 
and symmetrical, typically monovalent electrolytes were 
assumed (Mitchell and Ninham 1989; Winterhalter and 
Helfrich, 1989; 1992). 

Our calculations are limited to the equilibrium radius 
of a sphere or tube (with ends neglected), and we ig- 
nore the more subtle problem of determining the overall 
shape of a homogeneous vesicle with given volume, area, 
and leaflet area difference (a surrogate for the sponta- 
neous curvature) (Seifert and Lipowsky 1995; Mui 1995). 
Biological systems are very inhomogeneous; present are 
channels that span the membrane, and enzymes that cir- 
cumvent kinetic barriers to lipid repartitioning that are 
a feature of in vitro systems. In short, the electrostatic 
effects we are considering may be the physical adjunct to 
some of the biochemically defined actors, (wedge shaped 
lipids, adaptins, coatamers etc.), involved in vesicle and 
tubule formation. 

It should not be forgotten that electrostatic effects 
can be large in a biological context (Honig and Nicholls, 
1995); for a surface charge of -0.2|e|/nm 2 (corresponding 
to 10% of lipids each of size ~ 50A 2 possessing an elec- 
tronic charge), a screening length of lnm, and an aque- 
ous dielectric constant e w ~ 80, the surface potential 
eip ~ &bT, the thermal energy. Several parameters that 
occur repeatedly in this paper are diagrammed in Figure 
I for a single tubule of radius R on an isolated vesicle. 




FIG. 1. Schematic of a tubular section of a vesicle. The 
structure is assumed cylindrical and has radius R and thick- 
ness d. The interior and exterior screening lengths are k^ 1 
and Ko 1 , respectively. Our calculations assume L R. 

In the following section, we recapitulate the nonlinear 
electrostatic free energy (there is no unanimity among 
prior papers) and collect several formulae for the lin- 
earized free energy, and the leading 1/R term in the 
nonlinear free energy. We emphasize the effects of multi- 
valent electrolytes and interleaflet coupling. The surface 
charge responds to the curvature under two possible sce- 
narios; either not at all i.e., fixed charge per midplanc 
area; or maximally, with charge fixed per area of each 
water-lipid interface (and separate reservoirs assumed for 
each leaflet). Equilibrium radii are very different in the 
two cases. 

We also take literally the isolation between the inte- 
rior and exterior electrolytes effected by the membrane 
of a giant lipid vesicle in vitro to motivate a calcula- 
tion within a fixed number ensemble for the interior. 
Theories which assume a thermodynamic reservoir {e.g. 
Poisson-Boltzmann) are incorrect unless the total inter- 
nal charge is set to a particular value to within 0(ei/e w ) 
times the total surface charge. Deviations from this value 
can qualitatively change the conclusions. In the Results, 
we compute various quantities using realistic parameters 
and contrast the results and various assumptions. The 
concluding section reviews pertinent experiments with a 
view towards extracting values of parameters and sup- 
porting the various assumptions we have made. 

ELECTROSTATICS 

In this study, we assume all structures are smooth on 
atomic length scales and adopt the continuum limit. The 
problem is that of a cylindrical or spherical shell of thick- 
ness b — a = d = 4nm and dielectric constant eg — 2 
embedded in an ionic solution with dielectric constant 
e w = 80. We assume that changing the aqueous buffer 
ionic strength does not affect e w . These parameters ap- 
proximate a lipid bilayer membrane under in vitro and 
in vivo conditions. 

The total charge density p(r) of the system can be 
decomposed into an aqueous mobile contribution and a 
fixed piece, 

P(f) = X! PiJixed + Pi,mobile, (1) 
i i 

where the fixed charge contribution resides exclusively 
at the lipid-water interface and is characterized by a sur- 
face charge a. The surface charge results from ioniza- 
tion of e.g. phosphate, amine, and carboxyl groups on 
hydrophilic moieties of the lipids; these effective surface 
charges are typically negative. We do not consider spa- 
tially varying surface charges and lateral lipid phase sep- 
aration; however, for dilute charge densities, and smaller 
screening lengths, the interactions among charged lipids 
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are expected to be minimal (Hui, 1996) and the charged 
component would be approximately a uniform 2D ideal 
gas over the entire vesicle. Details of how a is affected by 
solution conditions, local membrane curvature, etc. are 
deferred until after results are presented. 

The mobile charge density is assumed to follow a Boltz- 
mann ensemble 



Pt,mobile{r) = Pi 



(2) 



where p°° is the bulk concentration of specie i and (3 1 = 
ksT. Substituting p mo biie[f) into Poisson's equation, 



V • (e(r)Vif) = -4Trp(r), 



(3) 



we obtain, for the electrostatic potential in the aqueous 
phase, 



£ 2 X7 2 ip = m sinhV> + n(z+)Y (e^ - 



e- z + 4 ') + 



(4) 



where efltp = ip, and m, n(z + ), and n(z—) are salt con- 
centrations in moles/liter of monovalent, z + -valent and 
z_-valent (e.g. for CaCl2, z + = 2; for Na2S04, z_ = 2) 
ions respectively. Lengths are expressed in nanometers 
and all numerical conversion factors are subsumed in 
I = 0.3083nm (T=300K). 

The electrostatic portion of the free energy for the lipid 
plus buffer system follows most readily in the limit where 
Eq. |] is linearized. Thus if a uniform surface charge is 
given, 



Gel = g a iVi 



(5) 



where i labels all surface charge positions. This prescrip- 
tion can be ambiguous in the nonlinear case (Sharp and 
Honig 1990). 

An alternative argument is then to use the expres- 
sion whose variation yields Eq. ^ and which agrees with 
the linear theory in the appropriate limit (Dresner 1963; 
Sharp and Honig 1990). Thus, with fixed surface charge, 



G e i = I dSaip - ^ 



<Tr 



U[<p] 



3L 

4-7T 



(6) 



d 3 r-\V<p\ 2 , 



where the primed integral is taken only over coordinates 
in the aqueous solution and the double-primed integral 
is taken over the bilayer region (a < r < b) occupied 
by the lipid acyl chains. This equation is also the large 
volume limit of a functional that was shown to be equiva- 
lent to the standard thermodynamic definition of G e i [a] 
(Dresner 1963), as well as expressions involving a pa- 
rameter integral (Marcus 1955). We have checked that 
Eq. H can be reduced to a form similar to the frequently 



employed charging integral (Marcus, 1955; Sharp and 
Honig, 1990), only in the one-dimensional case and for 
all dimensions in linearized theory. Furthermore, we have 
explicitly verified that the 1/R coefficient in an expan- 
sion of G e i(R) /2ttRL also agrees with that found from 
the charging integral when membrane leaflets are uncou- 
pled. 

The free energy with fixed surface potential is just Eq. 
with the first term omitted, G e e [ip(S)} — G e i [a] — 
J (p(S)adS. In all cases, Eq. || is understood to be eval- 
uated for a stationary solution, i.e, Eq. ^. 

The "potential" U [<p] is just 



{le(3) 2 U[ V } = m (cosh^- 1) 



+n{z_)~— ( e~* 





- 1 - 


1 


z+ 




H 


- 1 - 


1 


Z- 


Z- 



(7) 



in the units of Eq. |] and vanishes as ip — ► to make Eq. 
^ scale with the surface area. The surface value of ip is 
also to be varied and the net coefficient of 6<p(S) is the 
usual boundary condition 



<7 + -r-d n (p w - —d n ip e = 



(8) 



at all lipid- water interfaces. The normal derivative is 
taken positive from the lipid outwards. 



Linear Poisson-Boltzmann Solutions 

For potentials efif 1, the linearized form of the equa- 
tion of motion to be solved in the aqueous solution with 
the proper fixed charge boundary conditions is 



vV 



3 ^ = o 



(9) 



where Kj i0 are the effective inverse screening lengths in- 
side or outside the vesicle (z± > 0), 



£ 2 K 2 = m + ^n(z+)z + (z + + 1) + ^n(z-)z-(z- + 1). 

(10) 

Inside the bilayer (a < r < b), there are no charges, 
[7 = 0, and Laplace's equation holds. The cylindrically 
symmetric solutions in all regions are also explicitly dis- 
played in Appendix A. The corresponding linearized form 
of the free energy (Eq. ||) is G e £ = h J dSap and the free 
energy per midplane area, 2ttLR is 



1 



9ei = -z a aV>{a) 1 



d 

2R 



1 



a b (p(b) 1 



d 

2R 



(11) 



Expressions in the decoupled limit eg/(e w K,d) = have 
previously been found (Winterhalter and Helfrisch 1988). 
We will analyze consequences of the expansion 
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■<-'= c «+f +§+° (» 



(12) 



for screening lengths appropriate to physiological condi- 
tions. 

Further expanding the coefficients Cj in powers of 
e</(e w /«i), we obtain 



Co 









K dJ 



e e (Kidab - K da a ) 2 



(nidK dy 



O 



nd 2 



M 2 



(n d - 1) 



Cid - 1) 



e w K 3 tz 3 d 3 



M) 2 

0"a(K-K o ~ KiK 2 a - 2kI) + 



a a a b KiK (n - k») + 2ct & k + O 



6 ) m /S^tZ 



(13) 



(14) 



"neutral surface" (Petrov and Bivas 1984). We consider 
here the extreme case of conserved charge per midplane 
area, A = 2irRL. The physical reasons for this choice as 
opposed to fixed charge per area of membrane/solution 
interface is deferred to the Discussion. The only change 
to Eqs. 13, 14, and [IB] is a redefinition of a a ^- In the 
decoupled limit, upon replacing a a ,b — * &a,b/0- ± d/2R), 
the coefficients C, for g e e(R) become 



and 

C 2 



Trd 2 

< If 



ircP 



(Kid)'- 



:{l + Kid) 



(n d) 



;(l + n d) 



(16) 



M) 3 
3 



M) 2 

4 
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Kid 



(n d) 3 (K d) 2 
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K Q d 
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(17) 



with Co remaining unchanged. Note in Eqs. [l6] and [Tt] 

a a ,b are now the fixed charges per midplane area. Al- 
though C2 changes only slightly in magnitude, the be- 
havior of C\ is qualitatively different from Eq. 14. 



and 



C 2 



37rd 3 



( Ki d) s M 3 



O 



(15) 



With this notation, the total membrane bending stiff- 
ness is 2C2 + k m where k m represents bending stiffness 
from nonelectrostatic (such as mechanical) contributions. 
When C*i < the membrane spontaneously curves in the 
sense we have assumed (b = out, a = in) and it is interest- 
ing to understanding the physical origin of the effect. For 
nd <C 1 either <r& > a a (fy ~ K ) or k^ 1 > n^ 1 (<Jb — Co) 
favors bending. This prediction follows by replacing 
each screening layer by cylindrical capacitors; the larger 
charge, or thicker layer will prefer the exterior, i.e. the 
screening charge cloud expands. In the opposite limit 
(nd 3> 1), we can think of the two thin screening layers 
as flat; here, the energy is minimized by making the layer 
with the greatest energy/area be interior, i.e. occupy the 
side with less area. 

The influence of t£ is felt through the dimensionless 
combination ee/(e w nd) as expected from Eq. ^. This is 
always small under our conditions. Note that it is in- 
correct to estimate the energy as a triple of independent 
capacitors of thickness k" 1 ,^, and k~ x , i.e. o 2 d/et + 
0(a 2 1 Ke w ) which would make the lipid contribution ap- 
pear dominant. Also note that for ei/e w nd — > 0, C\ is 
antisymmetric with the interchange a a <-> Ub, K% <-> k - 
When coupling through the bilayer is not negligible, then 
the antisymmetry exists only for m = k or a a = 05. 

Finally, we consider a scenario in which the surface 
charges are fixed with respect not to their physical in- 
terface area, but with respect to an area defined by a 
radius in the interior of the lipid. This radius defines a 



Expansion of Nonlinear Poisson-Boltzmann 
Equation 

It is often noted (Mitchell and Ninham 1992; Winter- 
halter and Helfrich, 1992) that the first few terms of an 
expansion of G e i in terms of d/R or X/(kR) can be an- 
alytically obtained for the nonlinear theory. The first 
order (1/i?), is particularly simple to derive because one 
is expanding a stationary functional about its minimum 
and several terms cancel. The equation to be solved is 



vV 



dU 

dip 



dr 2 



r dr dip ' ^ 



with the boundary conditions given by Eq. Let cpo 
satisfy Eq. |l8| without the r~ 1 d r ip(r) term. Upon using 
the boundary condition tp — ► far from charged surfaces 
and U(0) = 0, this equation for po(r < a, r > b) can be 
integrated to yield 



2^o) 2 



U bo] = 0. 



(19) 



For negatively charged interfaces, po(r — a,b) is nega- 
tive and follows from Eq. [l9] with the boundary con- 
ditions [[e(r)V< / 9o]] = — 47rcr. Here, [[. . .]] denotes the 
discontinuity across the charged surface. The free en- 
ergy per midplane area from Eq. ^ is the sum g e g(R) = 
g+(R) + g~(R) + gi(R) where c/_, g+, and gi are contri- 
butions from the inner, outer and bilayer portions of the 
membrane system. One can readily show (Appendix B) 
that ge(R) contributes only even powers of 1/i?; the ef- 
fect of bilayer coupling via ee, however, is still felt through 
the surface values of ipo and d r ipo determined from the 
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boundary conditions. In Appendix B, we explicitly per- 
form the calculation of g e e(R) to 0(1/R) and obtain 



g ei ~ 3(00) + y (£_ - £+)z + 



d_ 

2R 



a b ip (b)- 



<r a ¥>o(a) 



47T 



/ (£_ - £+)dz 



(20) 



where all inner and outer "one-dimensional" energies 
£± = ^(d r ip) 2 + U [ip] are evaluated with ip = <po(z) and 
d r — > 9 Z . The second and third terms in Eq. give 
the 1/.R dependence of the free energy. We have also as- 
sumed that the inner radial surfaces of the tube (r = a) 
are far enough apart such that a» . 



Nonlinear Numerical Calculation 

For smaller R, higher order terms in 1/R must be taken 
into account or the full nonlinear solution for ip(r) and 
the corresponding free energy must be found numerically. 
The inner and outer solutions are matched via the sur- 
face charge dependent discontinuous derivatives on the 
lipid faces, R ± d/2. The two jump conditions furnish 
the two parameters in the solution to Laplaces equa- 
tion in the bilayer, ip(a < r < b) = Clnr + D. Inte- 
grating from r max - 5 > k , we introduce a parame- 
ter t such that ip(r max ) = tK (n o r max ) and ip'(r max ) = 
—tK Ki(K r max ) consistent with the linearized solution. 
The parameter t is then tuned until the solution at r = 
obeys ip'(0) — > 0. Solutions for ip(r) are then used in the 
nonlinear free energy Eq. ^. 
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FIG. 2. Electrostatic potential ip(r) (ksT/e) near nega- 
tively charged interfaces (a a = c& = — 0.2|e|/nm 2 ) centered 
at R = lOOnm with 50mM monovalent salt solution inside. 
The dotted and thin/thick solid curves are nonlinear Pois- 
son-Boltzmann theory for an exterior solution of 50mM mono- 
valent and 33.3mM divalent anions/cations, respectively. The 
long(short) dashed curves correspond to linear theory with 
50mM monovalent(33.3mM divalent) ions comprising the ex- 
terior solution. The linearised theory in the exterior region 
does not distinguish between cations (z+ = 2) and anions 
(*- = 2). 

Figure || plots the electrostatic potential ip(r) calcu- 
lated from both Eq. |J and its linearized form Eq. ||. 
For large radii (R 3> both the linear and nonlin- 

ear potentials ip(r) about the interfaces at R ± d/2 are 
essentially those of a one-dimensional theory and are in- 
sensitive to R. However, we will show that despite the 
similar behavior of tp, the behavior of g e i (R) for the linear 
and nonlinear cases are drastically different. 

The nonlinear free energies are calculated for various 
R and fitted for large R to polynomials in R . The 
1/R coefficients agree with those obtained by the analytic 
method of the previous subsection. The 1/R 2 coefficients 
are found to be drastically different from that expected 
from linear theory for all physiologically reasonable con- 
ditions as previously found (see Figure 1 of Winterhalter 
and Helfrich 1992). 



Net Charge/Finite Volume Effects 

Our calculations have assumed that the solutions on both 
sides of the bilayer are in chemical contact with large 
reservoirs. Applied to a closed vesicle, the interior po- 
tential is determined by the boundary conditions at the 
bilayer and the zero reference potential at infinity (out- 
side). Using Gauss's law on a surface inside the bilayer 
(a < r < b) we infer from the electric field calculated in 
the linear approximation a net interior charge (all solu- 
tion charges plus a a ) per area of 



€"W \^id Kod , 



0(e e /e w , 1/R)} 



(21) 



for the vesicle. Even though on a per volume basis this is 
very small, the consequences of imposing any other value 
inside, with the same order of magnitude, are dramatic. 
These concerns are not simply academic for once a vesicle 
is closed the interior charge may be fixed, so altering 
even the exterior electrolyte will lead to an imbalance 
between the charge present and that required under the 
open reservoir assumption. 

For a symmetric closed vesicle, the interior potential to 
within an additive constant can be found uniquely from 
the net charge in the solution; Gauss's law on the aqueous 
side of the bilayer (r = a_ ) gives the necessary boundary 
condition. To generalize Eq. ^ we have only to use the 
expression given by Dresner (1963) 
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U M = V N a in V- 1 P e-^ e ^d 3 r (22) 

£w a JO 

where q± — ±z±. By demanding stationarity under vari- 
ations in ip we find 



VV(r) = 



-Aire 



P _e/3u>\ 



_ e -e/3ifi(r) 



(23) 



where we have specialized to monovalent ions only. N± 
is the total number of cations/anions in the interior so- 
lution, n± = N±/V (with V the interior volume), (. . .) 
denotes V~ x J ■ ■ ■ d 3 r. Of course Eq. ^3] implies Gauss's 
law and is invariant under constant shifts in ip. Expand- 
ing the potential (p(r) — Sip+ (ip) about its interior aver- 
age value (ip) , determined largely by the exterior solution 
(ip(r — > oo) = 0) and the boundary conditions, we obtain 
W 2 ip(r) = K?0(r), 
where 



(p(r) = S(p 



(24) 



The solution to Sip(r) are similar to the potentials derived 
from Eq. ^ except for multiplicative factors and constant 
shifts, 



S(p(r < a) — 



2na (n + — nJ)e 
e w Ki h(Kia) 



Io(Kir) 



knT n + — ri- 



(25) 



from which we see (dtp) = 0. 

For purposes of illustration we evaluate G e i only in the 
case n + = n_ so the interior potential is constant and 
Sip = 0, a feature that persists in a nonlinear treatment. 
Here, (ip) ~ Aitdo a /ei + 0(ei/ e w nd) and is a factor of 
~ 40 larger than <p(S) calculated in the previous sections; 
however, in the linear limit with n+ = ri_, the free energy 
given by Eq. 11 still holds. The coefficients in a l/R 
expansion are 



Ci 



2nd 2 



O'aO'b — 



(26) 



and 



Co 



2nd 3 
Set 



Q / 9(cr a + a b f 



,(K d) 2 



8n d 



(27) 

Both C\ and Ci are a factor of ~ £ w /t£ larger than in 
the conventional theory. For small eg, the dominant term 
arises from a a tp(a) and its sign merely says that a charged 
shell prefers energetically to elongate. 



RESULTS 

In this section, we investigate the effects of charge asym- 
metry, ionic strength and multivalency of the inner and 
outer solutions on the linear and nonlinear electrostatic 
energies. A free energy of the form (Q) holds for each side 
of the impermeable membrane. Energies and charge den- 
sities will be expressed in units of feT/nm 2 and |e|/nm 2 
respectively. 
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FIG. 3. Linear and nonlinear free energies per area as a 
function of midplane radius for surface charges Gb = —0.15, 
and <Ja = —0.25. Dashed lines represent electrostatic contri- 
butions only, solid lines depict gTor(R) = g e i + 6k B T/R 2 . 
Linear solutions are shown in the upper pairs of curves while 
nonlinear Poisson-Boltzmann solutions are shown in the lower 
pairs of curves (nonlinear g e e(R) are fits to the open triangles), 
(a). 10~ 2,5 M ~ 3.16mM monovalent salt inside and outside 
the vesicle. For plotting purposes, the nonlinear energies have 
been shifted upwards by +0.8fc s T/nm 2 . (b). m = lOOmM 
on both sides. The nonlinear energies have been shifted by 
+0.01fe B r/nm 2 . 

The form of g e e{R) is displayed in Figs. || (a) and (b). 
Here we have for simplicity only considered the charge 
per physical leaflet area ensemble with a a = —0.25, <7b = 
—0.15, and monovalent ions both interior and exterior to 
the bilayer. In Fig. |(a), (m = 10~ 2 - 5 ~ 3.16mM), the 
linear free energy g e e(R) = g / cr{S)tp(S)dS (dashed line) 
is monotonically decreasing with increasing R, whereas 
the nonlinear g e i{R) has C\ < 0. When ion concentration 
is increased as in Fig ||(b), (m = lOOmM), nonlinear 
and linear results become similar in that C\ < 0. In 
fact, for the parameters used, the nonlinear g e i(R) shows 
no minimum in R implying that in the absence of other 
forces, the lipid tubule will collapse until the "hard wall" 
limit R ~ d is approached and the inner surface charge 
becomes a line. 

However, in addition to ionic forces, other shorter- 
ranged electrostatic and entropic interactions lead to me- 
chanical bending rigidities k m . The solid lines in both 
graphs plot g e i{R) + \k m /R 2 , with k m — 12/cgT. Un- 
der realistic physical conditions, the electrostatic contri- 
butions to the 1/R 2 terms in the free energy, C2 > 0, 
are small compared to measured mechanical k m val- 
ues, which in uncharged vesicles fall in the range k m ~ 
2 - mk B T (Song and Waugh, 1990; Andelman 1995). 
Therefore, the total free energy, for large R behaves as 

g TOT {R) ^C + ^+ (C2 + t" l/2) + O (1/i? 3 ) , (28) 

we will henceforth consider the 1/R 2 coefficient, C2 + 
k m /2 ~ k m /2 as an independently measured stiffness, of 
mostly nonionic contributions, and use an intermediate 
approximation, k m ~ 12k bT, where required. The bal- 
ancing of the 1/R and 1/R 2 contributions controls the 
size scales in bilayer bending and membrane structures, 
such that the free energy minimizing radius is 

(29) 

for C\ < 0, and the gain in free energy (relative to the 
flat state) at this radius is 

C 2 

gror (00) - groriR*) ^ ol , ^ onV ( 30 ) 

In Figs. S(a), (b), and (c), we plot the dependence 
of C\ on n\ for three different groups of surface charges; 



(cr a = -0.021, a h = -0.019), (cr Q = -0.21, a b = -0.19), 
and (a a = —0.25, 05 = —0.15), and both surface charge 
ensembles. 
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FIG. 4. C'i (fer/nm) as a function of monovalent 
ion concentration (logio scale, equal both sides). (a). 
G a = -0.021, G b = -0.019. (b). G a = -0.21, G b = -0.19. 
(c). G a = —0.25, Gb = —0.15. The filled circles and triangles 
indicate the approximations and parameters used in plotting 
g e e(R) in Figure [| Dashed(solid) curves show Ci under de- 
coupled(coupled) approximations. The upper set of curves in 
each graph (Ci > for large ni) correspond to surface charge 
per midplane area, while the lower branches (Ci < for large 
ni) correspond to surface charge conserved with respect to 
area physically occupied by charges. 
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In each figure, linear theory is represented by curves 
which sharply increase for low n\. The behavior of C\ 
calculated from nonlinear theory is drastically different 
for low ion concentrations, but approach as the linear the- 
ory at high concentrations. The pair of curves (C\ > at 
large n\) corresponding to conserved charge per midplane 
area remain positive. This is most easily seen in linear 
theory from Eq. |l6| and represents the dominant effect of 
increasing total interior surface charge. In the charge per 
midplane ensemble, the inner charge a a increases against 
a deeper potential (p(a) than o~b decreases; this electro- 
static work maintains C% > 0, thus biasing a negative 
radius of curvature (invaginations). Nonlinearities, as 
expected, mitigate these effects. Also shown is C\ in the 
decoupled limit (et/{e w Kd) — > 0), with dashed curves. 
We have checked that for all physically reasonable pa- 
rameters, g e i(R) varies by at most only a few percent 
when coupling between the two interfaces (at R — d/2 
and R + d/2) is neglected. 

The small effect of membrane coupling (for q = 
2, e w — 80) is clearly demonstrated in Fig. [| down to 
Tlx ~ O.lmM, especially at higher surface charges (com- 
pare Fig. |i|(a) with Figs. |J(b) an d (c)) because the differ- 
ence in surface potentials, ip(b) — <p(a), remains relatively 
small as the surface potential increase in magnitude non- 
lincarly. Since the coupling induced contribution of the 
bilayer enters through (p(b) — <p(a), its relative importance 
diminishes as \o~\ increases, particularly in nonlinear the- 
ory (see Appendix B, Eq. |3^). Also, note that since the 
interior and exterior solutions are identical, a symmetry 
exists, i.e. the interchange er Q «-> at leads to C\ — > — C\. 
The open triangles in Fig ^(c) indicate the parameters 
used in generating g e e(R) in Fig. ||. 

Figure H explicitly shows R* and <7tot (oo) — groriR*) 
based on Eqns. |29j and |0] with k m = l2k B T. Thus, 
a G\ of order O.lfcsT/nm is required to electrostati- 
cally induce radii of curvatures in the ~ 50nm range. 
For example, in the nonlinear cases plotted in Figs [|(b) 
and (c), at a 50mM monovalent salt concentration and 
conserved charge per leaflet area, Eq. 29 yields mini- 



on J(7nJ= 700mAf 



mum free energy tube radii of 182nm and 37nm respec- 
tively. The depths of these energy minima are given by 
Eq. H as 0.2xl0- 4 fc B T/nm 2 and 4.4 xlQ- 3 k B T /nm 2 
respectively. Similar values are obtained when charge 
per midplane area and the opposite charge asymmetries 
are considered. Coincidentally, the magnitudes of Ci are 
nearly equal for the two conserved charge ensembles when 
o = ipa + o&)/2 = —0.3; their magnitudes differ again as 
a > —0.3. Figure ^ shows that protuberances with radii 
relevant to biological systems can occur under appropri- 
ate conditions attainable experimentally in vitro and in 
vivo. 




0.000 
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FIG. 5. (a). Minimum energy radii R* (nm) (with 
the assumption k m = 12fc_eT) plotted as a function 
of charge asymmetry Aa = (cj — a a )/(a a + at) for 
mean charge a — (a a + <7(,)/2 = —0.1, —0.2, and 
-0.3|e|/nm 2 . Solid(dashed) curves represent fixed charge per 
leaftlet (midplane) areas. For fixed charge per midplane area, 
R* > corresponds to charge asymmetries opposite of those 
plotted, i.e., Act = (a a — ab)/(a a +at). (b). The associated 
free energy changes Ag(R*)(k B T/nm 2 ) (Eq. |^). Monova- 
lent salt concentration is fixed at lOOmM. 

The surface values, ip(S), where |v 5 (' , )| is maximal, are 
plotted in Figs. ^|(a) and (b) as functions of n\ and 
7i2, respectively in order to assess the validity of linear 
theory. Linear theory is expected to be accurate when 
zef3<p(r) < 1, although the linear C\ calculated in the 
charge per leaflet area ensemble seems to more accurate 
over a larger range of ionic strength than that calculated 
in the charge per midplane area ensemble (see Fig. ||). 
The effects of divalency on tp(S) are revealed in Fig ^(b). 
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Effects of divalent anions are quantitatively similar to 
those of monovalent salts, especially if we compare at the 
same concentration of monovalent cations which do most 
of the screening. However, for divalent cations, ^(S 1 )! 
is substantially reduced and linear theory is valid over 
a wider range of concentrations. The z + = 2 species 
is more effective at screening because they will balance 
more negative surface charge e/3ip(S) for the same penalty 
in entropy of mixing incurred. 




log[nj] 
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FIG. 6. Surface electrostatic potential at a negatively 
charged plane interface. (a). <p(S) as a function of 
monovalent salt concentration. Long dashed, dashed, dot- 
ted, solid, and thick solid lines denote surface charges of 
-0.02, -0.05, -0.1, -0.2, and -0.3|e|/nm 2 respectively, (b). 
Surface potential in the presence of a divalent salt solution. 
Thick(thin) lines indicate divalent cation(anion) solutions, 
and dashed(solid) lines correspond to a — —0.1 (a = —0.2). 



relative ionic strengths of the solutions, one can induce 
different screening and hence bilayer bending. For sim- 
plicity, we consider the extreme case of pure monova- 
lent and pure divalent salt solutions. Bilayer bend- 
ing due to solution asymmetry is demonstrated in Fig- 
ures 0(a) and (b), where C\ is shown as a function of 
a = —cr a = —o\, for both surface charge ensembles 
(charge/leaflet area, thin lines; charge/midplane area, 
thick lines) for interior monovalent ion concentrations of 
1, 10, and lOOmM from top to bottom within each triplet 
of curves. The necessary concentration of multivalent 
ions to achieve flaccid vesicles is assumed in the exterior 
solution, e.g. 2n^ ut ^ = [z± + l)n^ m ^(z±). Here, for di- 

valents, n^(2) = 2n[ out) /3 = .666, 6.66, and 66.66mM 
respectively. Figure shows that divalent cations (a) , are 
more effective at inducing larger curvatures (larger C\) 
than divalent anions (b). 
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Differences in the inner and outer buffer solutions can 
also affect membrane bending. By simply changing the 







FIG. 7. C\ {ksT /nm) as a function of a = — a a = — erf, for 
various solution asymmetries. Monovalent salt concentration 
in the interior is maintained at 1, 10, and lOOmM (dashed, 
dotted, and solid lines respectively) while an equiosmolar con- 
centration (0.66, 6.66, and 66.6mM) of divalents make up 
the exterior solution. Thick lines represent conserved charge 
per midplane area ensemble, (a). Exterior divalent cations, 
z+ = 2, (e.g. CaCb). (b). Exterior divalent anions, z_ = 2, 
(e.g.Na 2 S0 4 ). 

For small surface charges in the charge per leaflet area 
ensemble, linear theory, (Eqns. [l(] and |l4|) is expected to 
hold and yields positive(negative) slope for Ci(a) when 
ni is less(greater) than (3/2 + V2)(£/d) 2 M ~ 17.3mM 
However, at larger surface charges the behavior of C\ 
crosses over to a negative slope for this ensemble, espe- 
cially for higher salt concentrations (see Fig. |?j). Here, 
the nonlinear screening in the membrane exterior en- 
hances the decrease in | (p(b) | relative to that of \<p(a)\ 
as er is increased, thus increasing the electrostatic energy 
of the interior charge layer relative to that of the exterior, 
resulting in bending with C\ < 0. 



[n 1 ](in)=50mM 
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FIG. 8. Ci (fcsT jnm) as a function of n(z± — 2) in the 
exterior showing how substitution of divalent for monova- 
lent salt induces curvature. The interior solution is fixed 
at 50mM monovalent salt and the solutions are equiosmolar. 
Thin(thick) lines correspond to o a — <7b — & = — 0.1(— 0.2). 
The four lower curves (C\ < 0) correspond to fixed charge per 
leaflet area with solid(dashed) lines corresponding to divalent 
cations(anions). The four upper curves correspond to fixed 
charge per midplane area with other conventions unchanged. 

The values of C\ corresponding to conserved charge 
per midplane area (thick curves) increase quadratically 
with a agreeing with linear theory, until cr becomes large 



enough that nonlinear effects become important and sat- 
urates C\. The nonlinear effects are more prevalent for 
anions (Fig. @(b)) for the reasons discussed in relation 
to Fig. § 

Thus, we see that tuning the relative screening 
lengths is an effective way of inducing membrane cur- 
vature. Figure || shows C\ for surfaces of a a = <Jb — 
—0.1 (thin lines) and — 0.2(thick lines) as the screening 
in the exterior solution is scanned. Here, the vesicle 
interior is held at 50mM monovalent salt, while the 
exterior has a varying proportion of monovalent and 
divalent salt keeping the vesicle flaccid (this requires 
n (out) = o.050 -3n(° u *)(2)/2). The various divalent mix- 
tures required to induce positive curvature (Ci < 0, tube 
growth), or invaginations (C\ > 0) and the sensitivity of 
C\ to divalent concentrations are clearly shown. Since 
&a = exchange of the solutions on the two sides of 
the membrane interchanges C\ — > —C\. 
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FIG. 9. Ci (ksT /nm) as a function of a = — a a — —at 
for exterior cationic multivalencies, z+ =3(log dashed), 



4(dashed), 5(dotted), and 6(solid) lines, (a), 
(b). n{ out) = lOOmM. 



(out) 



= 10mM. 



10 



Finally, in Fig. ^|we plot C\ as a function of a = —a a = 
—Ob for vesicles with interior monovalent salt exterior 
higher multivalent ions (z± > 2). The influence of higher 
valencies on C\ is modest except at high surface charges 
or low total ionic strengths. 

DISCUSSIONS AND CONCLUSIONS 

In this paper we have presented calculations which sug- 
gest that electrostatic forces can control lipid membrane 
bending under realistic experimental and physiological 
conditions. Membrane deformations can be induced by 
an aqueous solution asymmetry between vesicle interior 
and exterior as well as by charge asymmetry between the 
two bilayer leaflets. Radii of curvature of the membrane 
bending can be in the neighborhood of 50-100nm typi- 
cally seen in biological processes such as vesicle budding, 
tubulation of ER and Golgi bodies, and endo/exocytosis. 

Tribulation and growth of necks from vesicles requires 
the membrane to nucleate such protuberances. Elec- 
trostatic changes in free energy alone can yield quali- 
tatively reasonable conditions for the formation of tube- 
like structures from a flat bilayer membrane. For exam- 
ple, if Ag TO T{R* ^ 40nm) ~ 0.0005£ B T/nm 2 , the to- 
tal free energy decrease is roughly O.ll/csT/nm length. 
Therefore, a flat membrane is stable against tube fluctu- 
ations of height L < lOnm. This critical length will be 
slightly greater due to the additional bending energy cost 
at the tube base; however, from experimental electron 
microscopy images (Mui et at, 1996), the base can be 
approximated with a portion of a torus of cross-sectional 
radius ~ lOnm. Using a bending rigidity of k m ~ YlksT, 
and the fact that this toroidal section has both positive 
and negative curvatures, we find that the total mechan- 
ical bending rigidity is qualitatively small and does not 
affect the energetics appreciably such that once a fluc- 
tuation exceeds ~ lOnm, it will continue to extend and 
lower g ee . 

Magnitudes of membrane surface charges experimen- 
tally measured indicate that our canonical estimate of 
— 0.2|e|/nm 2 is a reasonable physiological value. Sur- 
face charges measured in plant vesicles using parti- 
cle electrophoresis and dye fluorescence range from 
-0.03|e|/nm 2 to -0.24|e|/nm 2 (Sack et at, 1983; Chow 
and Barber 1980). These are averaged charge densities; 
higher concentrations of charged lipid could be recruited 
to incipient buds (phase separation) if electrostatics are 
playing a role. The numerous detailed chemical mecha- 
nisms of lipid-solvent interactions have not been modeled. 
The variation of lipid pK a 's with solution ionic strengths, 
the nonelectrostatic binding of cations to membrane sur- 
faces, and the hydrogen bonding among lipid hcadgroups 
can all affect the effective surface charge and is discussed 
by Tocanne and Teissie (1990). 




(c) 

FIG. 10. Microscopic models governing the distribution of 
charge and the neutral surface, (a). Bulky heads sterically fix 
charge per inner and out leaflet areas, (b). Small heads and 
entropically interacting acyl chains distribute charges at the 
heads according to midplane area. (c). Charge per midplane 
area preserved owing to interdigitated acyl chains. 

The response of surface charges to bending is also cru- 
cial in determining the electrostatic contributions to the 
free energy. The two extreme cases examined correspond 
to charge distributions which are li pid tail or head con- 
trolled and are depicted in Figs. [jlj. For large head 
groups where lipid packing is governed by steric and 
electrostatic interactions among head groups, the sur- 
face charge is approximately held constant as the bilayer 
curvature is varied. Conversely, when the tails occupy a 
thermodynamic area larger than the heads, the surface 
charge is approximately fixed with respect to the area de- 
fined by the midplane of the bilayer. Furthermore, other 
special phases may be important in determining charge 
distribution. For example, if the solvent contains glyc- 
erol or alcohols, the acyl chains of PC bilayers become 
interdigitated (Gennis, 1989) as shown in Fig. |l^(c). If 
this coupling is tight, the surface charge on the leaflets 
would be determined by the midplane area. 

It is interesting to note that the modes of charge con- 
servation mentioned above can be correlated with a filling 
parameter defined by Israelachvilli (1979), 

/ = v/ai, (31) 

where v is the effective volume of the lipid tail, a is 



11 



the thermal area of the head or the tail group near the 
head, and I is the vertical thermal length of the tail. 
For / ~ 1, the lipids are schematically represented by 
cylinders. For / > 1 and / < 1, the lipids behave mi- 
croscopically approximately as cones and inverted cones, 
respectively. The / > 1 ( "tail packed" ) or inverted cone 
/ < 1 ( "head packed" ) structures will probably be asso- 
ciated with fixed surface charge with respect to midplane 
area and fixed surface charge with respect to leaflet in- 
terfaces (r = a, b) ensembles, respectively. Depending 
on the chemical composition of the bilayer, an interme- 
diate neutral surface between R and a, b is also possible 
(Petrov and Bivas 1984, Gennis 1989). 

Correlating the cartoon in Figs. [l(] and the arguments 
above with bilayer chemistry may serve as an impor- 
tant guide in developing in vitro experiments and un- 
derstanding biological processes. Studies of the phases 
of concentrated lipid solutions suggest that / is related 
to lipid micellular, planar, and inverted micellular struc- 
tures. For example, lipids such as PC, PS, PI, PC, and 
Sphingomyelin at nonacidic conditions and in the ab- 
sence of divalent cations (Kates and Manson, 1984) form 
stable planar bilayers, implying that inter-lipid interac- 
tions through the tail and headgroups are comparable. 
Lysophospholipids, on the other hand, form micelles with 
the headgroups pointing outward into the aqueous phase, 
and can be modelled by inverted cones, with the cone 
apex at the midplane of the bilayer, commensurate with 
strong headgroup interactions. 

In general, factors that increase the effective acyl 
chain area relative to that of the headgroups favors 
H/7 (inverted hexagonal) phase over the L Q (planar) phase 
(Cullis et al. 1985). As in the case of surface charge, 
chemical conditions affect how lipids are packed in a bi- 
layer. For example, a low pH tends to increase headgroup 
association through hydrogen bonding, (Boggs 1984) and 
cations, such as Mg +2 and Ca +2 , also decrease head- 
group size by dehydrating them. This enhances the like- 
lihood of the H/7 phase. Lee, Taraschi and Janes (1993) 
have qualitatively applied the lipid shape concept to the 
POPtdEtn/PtdEth binary lipid mixture. Under their 
conditions, PtdEth has a large tail area which promotes 
the formation of inverted hexagonal phases in the lipid 
mixture (Lee, Taraschi, and Janes, 1993). In our study, 
vesicles containing similar types of lipids may be expected 
to conserve surface charge with respect to midplane area 
2ttRL. 

Probably the least accurately measured parameter in 
bilayer electro-mechanical models is the coefficient of 
1/i? 2 in gTOr{R)- Measurements using various methods 
have yielded scattered results. For example, in phospho- 
tidylcholines (PC), large values of k m are obtained when 
employing tube bending measurements, intermediate val- 
ues are obtained when a fluctuation mode analysis is per- 
formed, and low values are extracted when electric field 
deformation of spherical vesicles is analyzed (Andclman 
1995). For example, Song and Waugh (1993) measured 
k m ~ 28k B T for SOPC by mechanically pulling tethers. 



This value increased by ~ 3-fold when ~ 45% cholesterol 
was added. Mode fluctuation measurements by Mutz 
and Helfrich (1990) on lipid vesicles ~ 10/im yielded 
k m ~ 4k B T, 28k B T and l00k B T for galactosyldiglyc- 
eride, DMPC, and DMPC + 30% cholesterol bilayers re- 
spectively. Finally electric field deformation studies by 
Duwe et al. (1990) on DGDG and Egg Yolk PC vesicles 
yielded k m ~ 2k B T and k m ~ 5k B T respectively. All 
of these measurements were performed under neutral pH 
and low salt concentrations where electrostatic contribu- 
tion to the total bending rigidity fc m + 2C2 may be impor- 
tant. However, without careful experimental control of 
the solution ionic strengths (and surface charges), vari- 
ation in Ci can lead to the discrepancies reported. The 
enhancement of membrane bending stiffness upon surface 
association of uncharged polymers has also been studied 
by Evans (1996). Throughout this paper we have simply 
subsumed all these effects into an intermediate value of 
k m + 2C*2 — 12k B T, where C2 -C k m at the higher salt 
concentrations considered. 

Song and Waugh (1990) also measured bending stiff- 
nesses of artificial mixed POPS-SOPC vesicles as a func- 
tion of surface charge by varying the composition of 
charged POPS. They found no difference in total bending 
stiffness for , 2 , and 16% POPS (~ 35k B T). However, 
the solution ionic strengths and effective surface charges 
in these experiments were not precisely controlled and as- 
suming surface charges appropriate to the added amounts 
of POPS, even screening lengths ~ 30nm can cause C2 
to saturate at a maximal value <C k m (Winterhalter and 
Helfrich 1992). 

Experiments on artificial liposomes hitherto have not 
carefully considered the effects of solution ionic strength 
on deformation, tubulation and budding. In the experi- 
ments of Mui et al. (1996), a pH gradient across a bilayer 
was used to flip lipids with different pKa's across from 
one leaflet to the other. Tubulation was induced which 
the authors attributed to a leaflet area imbalance driven 
by the lipid exchange. However, as shown by Hope et 
al. (1989), pH induced transport of lipids across the bi- 
layer also changes the relative charges between the bilayer 
leaflets. These experiments therefore do not isolate the 
electrostatic effects we have calculated, even though, as 
we have shown, the electrostatic component may play a 
significant role in vesicle shape changes. Using natural 
Golgi bodies, Cluett et al. (1993) performed experiments 
with Brefeldin A (BFA) and approximately 50mM ion 
concentration. Tubules up to 7/iin in length grew when 
BFA was added to prevent binding of coat proteins and 
budding. The tubes were of the same size, ~ 70nm, as 
typical budding vesicles, suggesting that a common con- 
trolling factor such as electrostatics is not unreasonable. 

Besides continuum electrostatics, there are numerous 
other chemical and biological effects which can alter the 
mechanical properties of a bilayer. In particular, ion 
binding and hydrogen bonding effects have not been con- 
sidered. However, experiments hitherto have not care- 
fully controlled parameters affecting even the electrostat- 
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ics: solution ionic strength and surface charge. We pro- 
pose that in vitro experiments on large vesicles be per- 
formed under flaccid conditions with ionic strength as 
well as pH carefully measured. Systematically varying 
surface charge may also be appropriate. 

Our calculations have shown how adding neutral 
dopants to a bilayer can indirectly induce morphological 
changes of electrostatic origin. For example, if enough 
cholesterol, a rigid molecule with a relatively small head 
group, is incorporated equally in the bilayer leaflets, upon 
bending, phospholipid charges will be conserved with re- 
spect to a neutral surface closer to the the midplane ra- 
dius R as indicated by Figs. [lO|b, c. Although k rn > 
is also expected to change in magnitude, the sign of C\ 
is very sensitive to how charge is conserved. With asym- 
metry in er Qj b (Fig. ^) or interior/exterior screening (Fig. 
^), the addition of cholesterol can change the sign of C\ 
and determine whether tubules grow outward or invagi- 
nate. A gradient in lipid composition occurs biologically, 
for example with cholesterol in the Golgi-ER membranes 
(Bretscher and Munro 1993). Such membranes are con- 
stantly tubulating, budding, and recycling and their local 
cholesterol content may determine the neutral surfaces 
which govern the electrostatic component of these pro- 
cesses. 

We have also shown that provided a pure lipid bilayer is 
impermeable to ions, the conservation of charge in the in- 
terior of a closed vesicle can completely alter the electro- 
static energies calculated assuming thermodynamic reser- 
voirs. Only for a certain interior charge density specified 
to 0(a a /R) will the conventional result apply, though for 
a few tubes growing from a sphere of much larger area 
our charge imbalance analysis should be applied to the 
sphere only. For in vitro studies of pure bilayer vesicles 
it is experimentally difficult to control the interior charge 
to the required accuracy and thus large electrostatic ef- 
fects are likely to be present at least initially before the 
inner and outer solutions have equilibrated. 

For an interior charge of order <j a ,b per area the internal 
potential will be of order (e w / 'e^ksT '/ 'e, which is absurdly 
large for biological membranes, and even for an artificial 
bilayer will cause ions to traverse it. If clcctrogenic ion 
pumps act to maintain a electronic potential difference of 
several ksT across an organelle's membrane there may 
still be an effect on the optimal curvature comparable 
to what we have calculated in the Results section. Our 
calculations are most apt for a membrane where pores 
allow small inorganic ions to equilibrate while multivalent 
proteins are localized to one side and control the bending. 

Finally, we have verified that the effects of multivalent 
species (z a > 2) in the surrounding buffer solution are 
similar to those of the divalent solutions which we have 
treated in more detail. This has implications for vesicle 
budding assisted by adaptin/clathrin or dynamin pro- 
teins, whose molecular charges can interact and screen 
those at the bilayer surfaces. 
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APPENDIX A 

The solutions to the linear equation in the geome- 
try of Fig. [I] under appropriate electrostatic boundary 
conditions are displayed: 



<p(r < a) = 4iro- a K baKi(K b) 
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APPENDIX B 

Write the solution to Eq. |l8| as an expansion in 1/R; 
<P = fo + (fii where tpo is the solution to the one dimen- 
sional Poisson-Boltzmann equation for a flat interface, 
with boundary conditions appropriate to surface charges 
a a .b and cpi = 0(1/ R). Expanding G e t in for example 
the bilayer (a < r < b) and outer regions (r > b), 

%±^- ~ a b b(Mb) + Mb)) -^f Q(3^o) 2 + 
(d r ip )(d r ipi) + U [ip ] + tpiU' [<p ] )rdr- 



1L 

4-7T 



-(drtpoF + (0r¥>o)(0r¥>l) f *T + 0(1/ R) 



(34) 



Using the identity (d r (po)(d r ipi) = V- (ipi V(^o) — ^iV 2 ^, 
and recalling that to the order in R to which we are work- 
ing, we can approximate V 2 ^o — df,ipo = U' [ipo] in the 
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first integral, one sees that the U' terms cancel. The to- 
tal derivatives in both integrals reduce to surface terms 
which cancel <Jbbipi(b) and the corresponding term at the 
inner surface when G_ is included. This is no accident 
since the electrostatic boundary conditions just express 
the stationarity of (||) with respect to variation in the 
surface value of ip. In fact since d r tpo satisfies boundary 
conditions based on er Qi {,, (pi(a,b) = 0. Our argument 
does establish that tpo could satisfy the R — > oo bound- 
ary condition for the fixed charge per midplane area en- 
semble without affecting the 1/R coefficient. Of course 
the explicit cr a ,f, that occurs in G± should be correct for 
the particular ensemble. 

To extract the R dependence from the integrals over 
ipo we write rdr — > ±(i?± d/2 ± z)dz in G± and expand. 
Incidentally this argument shows that for a sphere the 
coefficient of 1/R is doubled, due entirely to the varia- 
tion of surface area with radius. The integral within the 
bilayer is symmetric about R and does not contribute to 
the 1/R coefficient. For numerical purposes it is conve- 
nient to use Eq. [l9] and express the remaining integrals 
in Eq. |2^ using <p , which is a monotonic function of r, 
as the independent variable, i.e., 
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